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Abstract 

We compute the proton-proton multiplicity distributions at LHC energies in the framework of a multiple 
scattering model assuming a Poisson distribution for each inelastic collision. Multiple scattering is essential 
to broaden the multiplicity distribution. We obtain approximate KNO scaling for small pseudo-rapidity 
intervals (I?;! < 0.5) and sizable KNO scaling violations for larger ones, in agreement with experiment. 



1 Introduction 

Multiplicity distributions in proton-proton collisions have been measured by the CMS collaboration at 
^/s = 900 GeV, 2.36 TeV and 7 TeV for central pseudo-rapidity intervals |r/| < 770 with r/o = 0.5, 1.5 and 2.4 
[T]. The measured rapidity distributions are much broader than Poisson distributions and show approximate 
Koba-Nielsen-Olssen (KNO) scaling [2] for the smallest pseudo-rapidity interval 770 — 0.5 with increasingly 
larger scaling violations as the length of the interval increases. For 770 — 2.4 the scaling violation takes place for 
z = n/ (n) > 3. Here n is the event charged multiplicity and (n) its average value in the considered 77 interval. 

Similar features have been observed at lower energies, between SPS and pp collider [3]. Actually, in this 
energy range the scaling violations for rjo — 2.5 start earlier {z > 2) and are numerically larger than in the LHC 
range. These features were well described in the framework of the Dual parton model (DPM) [3] and Quark 
gluon string model (QGSM) [5]. DPM and QGS are multiple-scattering models in which each individual inelastic 
collision is the superposition of two strings and the weights of the various multiple-scattering contributions are 
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given by a quasi-eikonal model (or perturbative reggeon field theory). One assumes a Poisson distribution for 
each string for fixed values of the string ends. The broadening of the distribution is due both to the fluctuation 
in the number of strings and to the fluctuation of the string ends. For energies in the LHC range and pseudo- 
rapidity intervals of limited length as the ones discussed here, the effect of the fluctuations of the string ends is 
negligibly small and DPM reduces to an ordinary multiple scattering model with identical multiplicities in each 
individual scattering |6j . The broadening of the multiplicity distributions is then entirely due to the fluctuation 
in the number of inelastic collisions. 



2 The Model 



Let us consider the charged particle per pseudo-rapidity unit dNPP/drj = (da^P / di]) / a^jj where the nu- 
merator is the charged single particle inclusive distribution and cr^^ the non-diffractive pp cross-section. At 
mid-rapidities and high energies, with identical multiplicities in each individual scattering, we have 

dNPP 



dNPP 
drj 



1 daPP 



'ND 



drj 



k>l 



k>l ' 



(k) 



dr] 



(1) 



Here (fc) is the average number of inelastic collisions and dN^^ /drj the charged multiplicity (per pseudo-rapidity 
unit) in an individual collision. Note that the general formula 54 for dNPP/ drj in DPM reduces to the simple 
expression ([T]) when all string contributions (and therefore all individual scatterings) are identical. This turns 
out to be the case at LHC energies for the small central rapidity intervals under consideration. dN^^ /dr] 
is then equal to twice the string multiplicity per pseudo-rapidity unit. This quantity is independent of s at 
mid-rapidities and high energies. 

Let us turn now to the weights 1^ for the occurrence of k inelastic collisions and their energy dependence. 
Following [6l [7] we use a quasi-eikonal model with exponential residues in t in which diffractive contributions 
are included as intermediate states in the eikonal model . We have [51 [7] 
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Here £, = £n{s/so) with so = 1 GeV^, crp = 8njp exp(AO and Z = 2Ce1p exp(A^)/(i?2 + a'p^). 

In (d) (jp is the Born term given by Pomeron exchange with intercept ap(0) = 1 + A. The (non-diffractive) 
inelastic cross-section in eq. (1) is cr^piiO = E c^feCf)- As for the numerator of eq. (f ) one obtains using (2) 
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This is a well known identity known under the name of AGK cancellation [S]. It implies that all multiple- 
scattering contributions vanish identically in the single particle inclusive distribution da /drj. Only the Born 
term contribution crp is left. 

The mid-rapidity values of dNPP /drj measured by the ALICE |9j and CMS [TO] collaborations have an s- 
dependence dNPP /drj{rj* = 0) ^ s°-^^. With cr^^ behaving approximately as s° °® [TT] we find from eqs. ([l]) 
to ([3]), Q;p(0) — I -I- A with A ~ 0.19 [6j, a value substantially larger than the one usually considered. In 



the following we take A = 0.19. The values of the other parameters in eq. ^ are [5] : a'p — 0.25 GeV^^, 
— 3.3 GeV^^, 7p = 0.85 GeV^^ and Ce — 1.8. The parameters and a'p control the i-dependence of the 
elastic peak and Ce contain the contribution of diffractive intermediate states. Its value is obtained from the 
ALICE results for single and double diffractive cross-sections :11 . The total pp cross-section is given by 

<Jtot{s)^ap f{z/2) , /(z)=^L^ (4) 

Note that eq. Q is obtained neglecting the real part of the scattering amplitude. However, only its imaginary 
part contributes to the weights {k > 1) and to (Tjvd -which are needed for the calculation of the multiplicity 
distribution. 

The charged multiplicity per unit pseudo-rapidity is given by 

dri dr] }J ak tr^''^ 
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As shown in [B] the model reproduces the energy dependence of dN^P/drj as well as its absolute values with 
dN^^ /d7]{ri* = 0) = 1.5 -consistent with twice the string multiplicity value. 

The results for ct^^ and dN^P /dr]{r]* — 0) at various energies are given in ref. [6]. Since they are needed in 
the calculation of the multiplicity distribution we have listed them, together with crf^^, in Table 1. 

Table 1: Values of the total and non-diffractive pp cross-sections and the charged particle pseudo-rapidity 
densities in the central rapidity region for the energy range between 200 and 50000 GeV. 



s (GeV) 


dNPP/driiy* = 0) 


f^ND (mb) 


<^Z (mb) 


200 


2.99 


31.22 


41.62 


540 


3.50 


38.97 


54.39 


900 


3.82 


43.33 


61.85 


1800 


4.34 


49.64 


72.93 


2760 


4.71 


53.77 


80.27 


5500 


5.42 


60.78 


92.92 


7000 


5.70 


63.33 


97.57 


14000 


6.61 


70.99 


111.57 


50000 


8.80 


86.19 


139.67 



3 Multiplicity distributions 



We turn next to the charged particle multiplicity distributions. The particle production in one-string process 
is expected to correspond to independent emission of clusters, described by a Poisson distribution. Since the 



convolution of Poisson distributions is also Poissonian, the multiplicity distribution for a process with k inelastic 
collisions will be given by a Poisson distribution -with average multiplicity equal to k times the one for a single 
collision, (or 2k times the one of an individual string) 

Here Uc is the number of emitted clusters and (n.c)o — (")o/^ average multiplicity in a single inelastic 
collision. The latter is equal to the average multiplicity of charged particles (n)^ divided by the average cluster 
multiplicity K, with 

d^^^2r;o^(r;*=0) = 3ry„ (7) 



valid at high energy for the relatively small values of rjQ under consideration 
The cluster multiplicity distribution is then given by 
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Eq. ^ is strictly valid only for infinitely narrow clusters. However it can be shown to be also true to a 
high degree of accuracy for clusters decaying according to a Poisson law [H]. A value of the average cluster 
multiplicity K ^ lA has been obtained in a Monte-Carlo simulation where the clusters are identified with a 
realistic mixture of directly produced particles and known resonances [13|. The same value K — lA allows to 
describe the charged multiplicity distributions in e+e~ and £p scattering [I^]. The same value of K is also needed 
in DPM to reproduce multiplicity distributions, long-range and short-range correlations in pp collisions up to 
i/s = 540 GeV 0]. As discussed in [T3], multiple scattering is responsible for long-range rapidity correlation 
which increase with y/s. 



4 Numerical results 

In the approach described above the pp multiplicity distributions can be computed with no extra parameters. 
It is convenient to plot them in the KNO form (i.e. (n) P„ versus z = n/ (n)). In this case the multiplicity 
distribution of final particles is identical to the one of clusters, i.e. 

^(z) = {n) P„ - (71,) P„^ (9) 

where z = n/ {n) = ric/ {uc). Here n is the number of emitted charged particles and (n) its average value 
(n) — (fc) (n)g ~ 3?7o {k). Also ric — n/K and (ric) — {n) /K. 

The results for ip{z) at two energies (-^/s = 900 GeV and 7 TeV) and two values of ?7o (770 = 0.5 and 2.4) 
are given in Fig. 1 and compared with CMS data Ij. The corresponding predictions at 14 TeV are also shown. 
We see that, in agreement with CMS data, the model has approximate KNO scaling between 900 GeV and 7 
TeV for the small pseudo-rapidity interval rjQ — 0.5 -up to a large value of z, (z ^ 6) where ip{z) has decreased 
by three orders of magnitude. For rjQ = 2.4, the KNO scaling violation starts much earlier (z ^ 2.5) and the 
multiplicity distributions in KNO form get broader with increasing energy. 
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Figure 1: The charged hadron multiphcity distributions in KNO form at y/s — 0.9 TeV (black sohd hne) and 
7 TeV (red dashed hne) TeV in two pseudorapidity intervals, \ri\ < 2.4 (left) and \ri\ < 0.5 (right) compared to 
experimental data [1]. Predictions at 14 TeV (blue pointed line) are also plotted. 

5 Conclusions 

We have computed pp multiplicity distributions at LHC in the framework of a multiple-scattering model 
(DPM). Multiple-scattering models do not obey KNO scaling. Indeed, the multiple scattering contributions, 
which give rise to long-range rapidity correlations, become increasingly important when s increases and since 
they contribute mostly to high multiplicities they lead to KNO multiplicity distributions that get broader as 
s increases. On the other hand, the Poisson distributions in the individual scatterings lead to short-range 
rapidity correlations and give rise to KNO multiplicity distributions that get narrower with incresing s. Due 
to the interplay of these two components the energy dependence of the KNO multiplicity distributions (or of 
its normalized moments) depends crucially on the size of the rapidity interval |16| . For large rapidity intervals 
the multiple-scattering effect dominates and KNO multiplicity distributions get broader with increasing s. For 
small intervals the effect of the short-range component increases leading to approximate KNO scaling, up to 
z ~ 6. We have shown that the above features are maintained up to the highest LHC energy and that for a 
given pseudo-rapidity interval (?7o — 2.4) the rise of the KNO tail starts at a value of z that increases with 
energy. 
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